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PACS 73 . 23 . -b - Electronic transport in mesoscopic systems 
PACS 73 . 2 1 . La - Quantum dots 

PACS 72 . 15 . Qm - Scattering mechanisms and Kondo effect 
PACS 73.23.Hk - Coulomb blockade; single-electron tunneling 

Abstract. - We show that the conductance of a quantum wire side-coupled to a quantum dot, 
with a gate potential favoring the formation of a dot magnetic moment, is a universal function of 
the temperature. Universality prevails even if the currents through the dot and the wire interfere. 
We apply this result to the experimental data of Sato et al. [Phys. Rev. Lett. 95, 066801 (2005)]. 



Nanodevices owe much of their development to the 
theory of many-body phenomena. Consider, e. g., the 
single electron transistor (SET), a quantum dot bridg- 
ing two otherwise independent two-dimensional electron 
gases [1,2]. The competition between the Coulomb block- 
ade, which bars transport through the dot, and the Kondo 
screening of the dot magnetic moment by the electron 
gases [3], which favors low-temperature conduction, was 
discussed on blackboards [4] a decade before it surfaced 
in the laboratory [1] . By the time the first device was de- 
veloped, quantitatively accurate theoretical results were 
available. Chiefly important was the universal conduc- 
tance curve Gg ET (T) for the symmetric Anderson model 
[5,6], which was shown to match the temperature depen- 
dence of the zero-bias conductances in SETs and analo- 
gous devices. 

More recently, experiment has leaped ahead of theory. 
The development of complex structures, such as the side- 
coupled device [7-12], has motivated only qualitative pre- 
dictions. As Fig. 1 shows, the current in the side-coupled 
device is carried by electrons that can either traverse the 
quantum wire or hop to a quantum dot to skip the central 
section [12]. A Fano parameter q [13], defined below, mea- 
sures the amplitude for the latter process relative to that 
for the former. The limit q — > oo emulates a SET. For 
smaller q's, the wire bypasses the Coulomb blockade and 
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Figure 1: Side-coupled device. The gate potential V g controls 
the energy Ed of the quantum-dot level Cd- The open circle 
depicts the Wannier orbital /o. 



allows high-temperature conduction. Below the Kondo 
temperature Tk, the screening of the magnetic moment 
enhances the electronic flux through the dot and allows 
interference with the flow along the central portion of the 
wire [14,15]. 

Attentive to the diversity of experimental findings, 
we have applied numerical renormalization-group (NRG) 
tools to an Anderson Hamiltonian modeling the side- 
coupled device. The resulting essentially exact numerical 
data for the temperature-dependent conductance G q {T) 
will be detailed elsewhere [16]. Here, we focus the relation 
between G q {T) and the universal curve Gg ET (T). Our 
central result covers the Kondo domain, the set of dot en- 
ergies and dot-wire couplings that favor the formation of 
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a Kondo screening cloud. 

Since the thermodynamical properties of the Kondo 
crossover are universal functions of the temperature 
scaled by Tk [2,17-25], a mathematical relation between 
G q (T/T K ) and the universal function Gg ET (T/T K ) is 
hardly surprising. Nonetheless, the diversity manifest 
in conductances that, in contrast with Gc ET (T), rise 
with temperature and in conductance profiles (fixed- T 
conductance-vs. gate-voltage plots) that display antires- 
onances [12] rules out a proportionality between G q (T) 
and Gg ET (T). Instead, we will show that a linear map- 
ping binds the two functions: 

G q {^) - y = (gIeAt^) - y) cos2<5, (1) 

where S is the ground-state phase shift of the wire elec- 
trons, and G2 = 2e 2 //i, the quantum conductance. 

While the mapping (1) is universal, the phase shift and 
Kondo temperature are model-parameter dependent. At 
fixed temperature, the Fano parameter q controls the func- 
tional dependence of the conductance on the gate voltage 
V g . As q grows, valleys in the conductance profiles evolve 
into plateaus, a result in qualitative agreement with mea- 
surements. Most importantly, Eq. (1) affords quantitative 
comparison with experiment; as an illustration, we will 
present curves that reproduce the temperature-dependent 
conductances reported by Sato et al. [8] ; show that the dot 
moment was fully screened; and extract Tk and 8 from 
the data. Our results justify mathematically the authors' 
phenomenological treatment of their results. 

Overview. Preliminary to the formalism, we present an 
overview of conduction in the side-coupled device. We con- 
sider weak coupling to the wire, so that the dot occupation 
rid is nearly conserved, and illustrate the discussion with 
NRG plots of the temperature-dependent conductance. 

The device has three characteristic energy scales, set 
by (i) the coupling to the wire, which broadens the dot 
levels; (ii) the electrostatic barrier A at between adjacent 
dot occupations rid = N — 1 and rid — N; and (hi) the 
Kondo temperature Tk, below which the wire electrons 
screen the dot moment. The first two scales catch the eye 
in conductance profiles. The third one defines the thermal 
regimes T ^> Tk and T -C Tk displayed schematically in 
the top and bottom panels of Fig. 2a, respectively. 
Left: q > 1. If T > T K (top), the Coulomb blockade 
impedes conduction; to defeat it, the gate voltage must be 
raised, so that the rid = N and the rid = N + 1 ground 
states are nearly degenerate. The conductance profile is 
hence a sequence of narrow resonances. Upon cooling 
(bottom), little changes if V g makes rid even. For odd rid, 
however, the Kondo hybridization between the wire and 
the dot states allows conduction. The conductance G(V g ) 
alternates between insulating valleys (G = for even rid) 
and Kondo plateaus (G = G2 for odd rid) [27]. 
Right: q = 0. The pattern is reversed. If T 3> Tk (top), 
except at the resonant voltages, the flux through the wire 
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Figure 2: Bird's eye view of conduction through the side- 
coupled device, (a) Condutance G as a function of applied gate 
voltage for three representative Fano parameters q at temper- 
atures high (top panels) or low in comparison with the Kondo 
temperature Tk- (b) NRG results for G(T). The crosses, 
squares, and circles in each panel represent conductances cal- 
culated at the voltages indicated by the same symbol in the 
panels directly above it. The solid lines depict Eq. (1), with 8 
extracted from the low-energy fixed-point eigenvalues [26]. 



is ballistic. On resonance, the strong coupling to the dot 
blocks conduction, and G(V g ) dips to zero [28]. For T <C 
Tk (bottom) and odd rid, the screening cloud blocks the 
wire. The conductance thus alternates between ballistic 
plateaus and Kondo valleys [14]. 

Center: q = 1. There are now two conduction paths. 
If T ^> Tk (top), Fano antiresonances near the resonant 
voltages signal interference between the two currents. Off 
resonance, the electrons flow only through the wire, and 
the conductance remains close to £2/2. For T -C Tk (bot- 
tom) and odd nd, again the Kondo hybridization to the 
wire allows conduction through the dot. The Kondo cloud 
nonetheless blocks conduction through the wire, so that 
the off-resonance conductance is again close to Q2/2. 

For each voltage identified by two crosses, squares, or 
circles in Fig. 2a, Fig. 2b displays our NRG results for the 
conductance as a function of T normalized by Tk — the 
temperature at which the conductance is half the quantum 
conductance, G q {T = Tk) = G2/2. The agreement with 
the solid lines representing Eq. (1) is very good: even at 
the limits of the Kondo regime (crosses and circles, T = 
50 Tk) the absolute deviations are smaller than 0.01 e 2 /h. 

For q ^> 1 (left, 5 ~ 0), the device mimics a SET, 
the mapping (1) reduces to an identity, and the conduc- 
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tance decays from ballistic to zero along the universal 
curve Gg ET (T). The opposite extreme, q = (right, 
S w 7r/2), reverses the pattern: ballistic conductance at 
T 3> T K , and perfect insulation at T = 0. As intuition 
would dictate, and as general [8] and mathematical [29] 
arguments suggest, the temperature dependence comple- 
ments the universal function: G q (T) = Q 2 — Gg ET (T). 
For q r~j 1 (center, 5 7r /4), three voltages V g are 
highlighted. In each case, the conductance interpolates 
monotonically the high- to the low-temperature limits in 
Fig. 2a. For 2ed + U = 0, in particular, Fig. 2a shows that 
G q (T > T K ) = G q (T = 0) = g 2 /2; the conductance must 
therefore be constant, G q (T) = £2/2, and this is what 
Eq. (1) predicts for 5 = 7r/4. 

Model. Hamiltonians describing side-coupled nanos- 
tructures, in the geometry of Fig. 1 [12,14,28-30], or other 
arrangements [31-33] have appeared in print. While our 
analysis could start from any of the former, to keep the 
presentation self-contained we define the alternative three- 
component Hamiltonian H = Hd + H w +H w d- Here, Hd = 
£dnd+Undiridi, models the quantum-dot, where rid = c\cd 
is the dot-level occupation; Ed is the dot energy, controlled 
by the potential V g ; and U is the Coulomb repulsion. The 

second term, H w = J^k e kc[ c k + ( K / N ) J2kk' c l Ck 'i mod " 
els the wire, of length L. The energies ek, measured from 
the Fermi level, form a symmetric, structureless, half-filled 
conduction band of width 2D comprising N levels sepa- 
rated by the splitting A = 2D/N [34]. The scattering K 
describes the (fixed) gate potential applied to the wire [12]. 

The last component, H wd = ^ fe (Vfc/ 'VN)c\c d + H. c, 
couples the wire to the dot. The relatively small thermal 
energies ksT <C D authorize retention of the two lead- 
ing terms in the expansion of the coupling Vk in powers of 
€k [35]: Vk ~ Vk F +€k(dV/dek)k F - We rewrite this approxi- 
mation as Vk = V+npqVek, where p = 1/2D, to define the 
coupling V and the Fano parameter q [13]. Following NRG 
tradition, we introduce the shortand / = J2 k Ck/VN. 
For future reference, we note that, in this notation, the 
(particle-hole) symmetric (q = K = 2e d + U = 0) Ander- 
son Hamiltonian reads 

#o 5 = €kC l Ck + V ^ Cd + H - c ") - T( nd T - n di ) 2 - (2) 
k 1 

The conductance is more easily computed on a ba- 
sis constituted of fo and N — 1 other conduction op- 
erators a p = ^2 k a p kCk, where p — 2irn p /L [n p = 
-N/2 + l,...,N/2 - 1], such that {$,a p } = 0, and 
that {a p ,a p i} — 5 PP > [36]. With the shorthand /1 = 
J2 P d p /\/N — 1, the model Hamiltonian becomes 

H A =Y,^4 a p + K fofo + ( t o.fld q + Vf^d + R. c.) + H d , 
p 

(3) 

where N q = y/l + (npVq) 2 ; N q d q = f + npqVcd] to = 
DAf q /V3: and the energies l k are the conduction energies 
e k phase shifted by ir/2, i. e., e k = e k — A/2. 



If K — ► 00, the scattering potential decouples fo from 
the other states, freezes its occupation at /<j/o = 0, and 
forces the current through the dot. The condition 2e d + 
U = reduces Ha to the symmetric- SET Hamiltonian 

H SET = J2 i P a l a P+ V i(fi c d+ B - c.)-^> dT -n di ) 2 , (4) 
v 

where V q = n q V/(2y/3). Given the analogous definitions 
fo ~ J2k Ck and h ~ J2 P a pi we see tnat OIU y the phase 
shifted energies Ik = £k — A/2 distinguish H§ ET from Hq . 

The first term within parentheses on the right-hand side 
of Eq. (3), which couples the conduction states /1 ~ a p 
tod q <~ fo+TpqVcd, is the mathematical expression of the 
two conduction paths in Fig. 1 , one of which runs through 
the central wire-orbital fo, and the other, with relative 
amplitude npqV, through the dot orbital c<j. To flow 
from the left to the right side of the wire, the current must 
traverse d q . Accordingly, the Linear- Response Theory [37] 
shows that p q , the spectral density for the operator d q , 
controls the zero-bias conductance: 



c q (T) = g 2 N q 2 

J-D 



Pq{t,T) 



9_m 

de 



de. (5) 



Here /(e) is the Fermi function, and 
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Z ™„ f( e ™n) 



\(m\dl\n)\ 2 6(e mn -e), (6) 



where Z is the partition function, and e mn = E m — E n . 

Analysis. We are chiefly interested in the Kondo regime, 
the set of temperatures and model parameters favoring 
unitary dot occupation, i. e., such that the energies \e d \ 
(ed + U) to remove the dot electron (add a second electron) 
dwarf the energy fc^T and dot-level width T = npV 2 . The 
Kondo temperature Tk then sets the energy scale. 

At high temperatures, T » Tk, the Hamiltonian (3) 
lies close to an unstable local-moment fixed point (LM), 
which comprises a noninteracting spin-1/2 variable (the 
dot spin) decoupled from a conduction Hamiltonian [26] 



(7) 



where the effective scattering potential Klm depends on 
the model parameters. For Ha — Hq {Hg ET ), in partic- 
ular, K LM = (K LM — > 00). 

Diagonalization of the quadratic form (7) yields N 
eigenoperators gk and eigenvalues ek — tk — SlmA/tt: 



h 1m ~^2 £ kglgk- 



(8) 



Near the Fermi level, the phase shifts Slm are uniform [26]. 

In analogy with the definitions fo ~ J2 k Ck and fx ~ 
J2k ekCk > we can define the mutually orthogonal combi- 
nations of eigenoperators <f>o — Sfefl'fe/V^V and (pi = 
y/3/N ^2 k (ek/ D)gk- For any pair of constants ozq and a.\, 



ao/o + ai/i = (3o4>o + Mi 



(9) 
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where 0o and j3\ are linear combinations of a and a\, 
with coefficients fixed by Slm- 

Tabic 1 collects results for the symmetric [Eq. (2)] 
and the symmetric- SET [Eq. (4)] Hamiltonians. While 
columns 3-5 follow from the definitions of g^, 4>o and <j>i, 
columns 5 and 6 require explanations. For H A = H§ (i. e., 
q = K = 2e<i + U = 0), the operator d q reduces to f = <fi . 
Its spectral density p q = p^ is therefore equal to p^ a . 



Table 1: Properties of the symmetric Hamiltonians 





0~LM 


9k 


00 




P4>i 

P4 "' Mq 2 


H, s 

H§ET 




tt/2 


Ck 


fo 

h 


Dfi 


Po 

PSET 



For Ha = Hg ET , d q reduces — irpq V Cd/M q . To relate 
its spectral density to p^, choose two eigenstatcs |m), \n) 
of H§ ET with E m , E n w ksT <C D. In the identity 

HE^|n) = ^ + V, W »>, (10) 

the left-hand side is then much smaller than each ma- 
trix element on the right. This shows that (m\(f>i\n) w 
— -Kpq V(m\cd\n), a result equivalent to (m\(f>i\n} = 
—N q {m\d q \n) , and hence to the last column in Table 1. 

As the system is cooled, the wire electrons screen the 
dot moment, and Ha crosses over from the LM to a stable 
frozen-level fixed point (FL). The latter is but the conduc- 
tion band resulting from letting Klm — * 1/ '(t 2 P 2 ^-lm) in 
Eq. (7). Diagonalization of the FL Hamiltonian yields TV 
eigenvalues = — SA/tt, where, in conformity with 
Friedel's sum rule [38], the phase shift S = 5lm — tt/2. 



Table 2: Conductance at the two fixed points 



Fixed 


Phase 






G%ET 


point 


shift 




(J = tt/2) 


(6 = 0) 


LM 


(5 + 7T/2 


G 2 sin 2 S 


Q2 





FL 


6 


Qi cos 2 S 








The fixed-point physical properties are independent of 
temperature and energy. In particular, the LM and FL 
conductances are trigonometric functions of the phase 
shift. The expressions, which result from an extension 
of Langreth's argument [38], are recorded in Table 2. 

The effective antifcrromagnctic interaction Hj = 
JS-(Tav4>QtJ,4>OL> [18,26,39], between the dot spin S and 
the spin of the localized orbital <j> drives the Hamiltonian 
from the LM to the FL. The NRG trajectory is universal: 
scaled by UbTk, the eigenvalues of Ha are universal, and 
so are the corresponding eigenstates on the basis of the 
{gk} [17,26,40]. The spectral densities p^, and p ( p 1 in Ta- 
ble 1 are therefore universal functions of the ratios T/Tk 



and c/UbTk ■ To highlight these findings, we define the 
scaled energy £ = e/ksTx and temperature T = T /Tk- 

Next, we turn to the asymmetric Kondo-domain Hamil- 
tonians. To relate the conductance G q to the universal 
function Gg ET , we add to H A an infinitesimal harmonic 
perturbation, frequency u>, coupling the spectrum of Ha 
to an auxiliary orbital d at the Fermi level: 

H v = rid i d q e- i " t + H. c. (11) 

The golden rule shows that the spectral density p q is 
the response function for the thermally averaged transition 
rate j d—>A induced by H n : 

(jd^))T = (W/n)p q (nw,T). (12) 

Consider, then, the perturbative effects of H v upon the 
Kondo crossover [18]. Close to a (Fermi-liquid) fixed point 
H* , that is, for H = H* + SH, one can always construct an 
effective Hamiltonian H e f f that reproduces the spectrum 
of H to linear order in SH. Here, to follow a pedestrian 
route, we subject the sum Ha + H v to the Schrieffcr- Wolff 
transformation: H\ !S + H^f = e s (H A + H v ) e - S [39]. 
Besides substituting the spin-spin interaction Hj for the 
dot- wire coupling H w d, this casts Eq. (11) in the form 

H e /f = (77/^)3 Wd + ao/o + ai/i)e- iwt + H. c, (13) 

where ad, an, and a\ depend on V, £<j, U, and K. 

The unperturbed effective Hamiltonian H A ^ commutes 
with rid. Of its eigenstates, only those with = 1 are 
energetically accessible at the LM. It is safe to disregard 
the perturbation proportional to ad in Eq. (13), which 
couples them to the subspaces rid = and rid = 2. 

The other two perturbations conserve rid and are impor- 
tant even at very low energies. It is convenient to project 
them upon (j> and <j)\ , because the spectral densities for / 
and /i are phase-shift dependent. Substitution of Eq. (9) 
for a f + ai/i brings Eq. (13) to the form 

Htf' = ( V /Af q )dUPo <Po + Pi 0i ) + H. c. (14) 

In the (LM to FL) crossover, just as the Kondo Hamil- 
tonian H^f = e s HAe~ s is equivalent to Ha and the 
diagonalization of H e J^ yields the physical properties for 
Ha, the effective perturbation H*ff is equivalent to H r] 
and application of the golden rule to H*ff yields the tran- 
sition rate induced by H v . Out of the terms then resulting 
from Eq. (14), only those proportional to Kml^Jln)! 2 and 
|(m|</4|n)| 2 contribute to G q [41]. The cross terms disre- 
garded, the last two columns in Table 1 lead to 

(j-^ A )T = ^(^p S 0+(3 2 p S SET)- (15) 

Comparison with Eq. (12) then shows that 

Pg (£,T) = (p 2 /JV q 2 )p s (£,T)+Plp s SET (£,T), (16) 
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and Eq. (5) yields an expression for the conductance: 

G q {T) = Pl G S (T) + M q 2 (3 2 G S SET (T). (17) 

Substitution of the expressions in the LM line of Table 2 
for G q , Gq, and Gg ET shows that 0$ = sin 2 8, while the 
expressions in the FL line show that N q {3\ = cos 2 8. With 
this, Eq. (17) becomes 

G q (T) = G$(T) sin 2 8 + G S SET (T) cos 2 6, (18) 

and to complete the derivation of Eq. (1) we only have to 
recall that Gq = Qi — G% ET . 



Table 3: Parameters for the five runs in Fig. 2b. U = 0.20 D. 



-e d /D 


V/D 


-K/D 


q 


S/ir 


k B T K /D 


0.10 


3.6xl0~ 4 


100.0 


100 


0.00 


1.3xl0~ 5 


0.17 


0.021 


0.315 


1 


0.22 


l.lxlO- 5 


0.10 


0.021 


0.315 


1 


0.25 


l.lxl0~ 9 


0.022 


0.021 


0.315 


1 


0.28 


l.lxlO" 5 


0.10 


0.056 


0.000 





0.5 


1.4xl0" 5 



Discussion. The phase shift 8 and Kondo temperature 
Tk in Eq. (1) depend on the model parameters. Almost 
invariably, they have to be computed numerically, because 
the perturbative expressions for 8 and Tk are accurate 
only in corners of the parametrical space [3,17,20,26,42]. 
Exceptions are the particle-hole symmetric Hamiltonians 
H$ [Eq. (2)] and H§ ET [Eq. (4)], for which 8 = vr/2 and 
8 = 0, so that Eq. (1) reduces to G q (T) = Q 2 -G S SET (T) = 
Gq(T) and G q (T) = G§ ET (T), depicted in the right and 
left panels of Fig. 2b, respectively. 

In other regions of the Kondo domain, 8 lies between 
— tt/2 and ir/2. In the Kondo crossover, the conductance 
on the left-hand side of Eq. (1) changes by less than Q 2 , 
from Q 2 sin 2 8 to G2 cos 2 8. The rise or decay of G q (T) is 
centered at Q2/2 and proportional to G§ ET — Qij2. 

The central panel of Fig 2b shows examples. The solid 
lines are Eq. (1) with 8 (Table 3) extracted from the FL 
single-particle eigenvalues in the NRG diagonalization of 
the pertinent Hamiltonian Ha] and Tk (Table 3), from a 
fit of the universal magnetic susceptibility [18, 20] to the 
susceptibility computed for the same Hamiltonian. No 
adjustable parameter is therefore involved in the excellent 
agreement between Eq. (1) and the NRG data for the con- 
ductance. We have applied the same procedure to more 
than 100 NRG runs sampling the Kondo domain; in each 
one, the agreement was equally good. 

In brief, the exact universal mapping (1) interpolates 
between G{T > T K ) and G(T < T K ), While the two 
extremes, described by single-particle Hamiltonians, are 
accessible to a variety of techniques — e. g., the Landauer- 
Buttiker formula or scattering-matrix analyses [43] — , the 
interpolation covers the temperature range beyond the 
reach of simple analyses. 




100 200 300 400 500 600 700 
T(mK) 

Figure 3: Comparison with experiment [8]. The crosses and 
squares are the conductances G at the indicated gate potentials 
V g ; the temperatures measured below lOOmK carry large un- 
certainties [8]. The dash-dotted and solid lines depict Eq. (19). 
The inset shows that, for each V g , the appropriate Tk straight- 
ens the plot of G vs. G s and yields Eq. (19). 

Comparison with experiment. To allow for the in- 
evitable background currents, we exploit the linearity in 
Eq. (1). Given a set of N experimental pairs {Gi, T;}, from 
a trial Kondo temperature T K we generate the dimension- 
less temperatures % — Ti/T K (i — 1,...,N). Next, we 
invert the function G s (T) to determine the universal con- 
ductance G si for each %. If the plot of Gi vs. G si is 
straight, T K is the Kondo temperature. If it is not, we it- 
erate. Visual inspection (numerical evaluation of the cur- 
vature) determines Tk within 10% error (to an accuracy 
limited only by the experimental dispersion). 

As an example, the inset of Fig. 3 treats the conduc- 
tances in Fig. 3b of Ref. [8] , measured with the gate po- 
tentials V g = — 795 mV (squares) and — 775 mV (crosses), 
and yields the Kondo temperatures Tk = 850 mK and 
1230 mK, as well as the linear relation 

G(T) = l.9e 2 /h - Q.5G S SET (T). (19) 

For T > T K , G s set (T) -> 0, so that G -> 1.9e 2 /h, close 
to the measured off-resonance conductance (1.8e 2 /h) [8]. 

The agreement with the solid and the dash-dotted lines 
representing Eq. (19) in Fig. 3 shows that the data are 
in the Kondo regime and the contact resistance is neg- 
ligible. There is, however, a background current: at 
T = 1, Eq. (19) yields G(T = 1) = \Ae 2 /h; the ex- 
cess Gfc = 0Ae 2 /h over G(T) = e 2 /h, predicted by 
Eq. (1) is the background conductance [8]. If we define 
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G = G — Gb, then Eq. (19) takes the universal form 
G - e 2 /h = —0.5(G§ ET - e 2 /h), from which we find the 
shift: cos 26 = -0.5 (6 w tt/3). 

In conclusion, we have shown that, in the Kondo regime, 
measured from C? 2 /2 and scaled by cos 26, the conductance 
of side-coupled devices is a universal function of T/Tk- 
Our application to experimental data identified inequiv- 
ocally the measured thermal dependences with Kondo 
screening; detected a background current; and determined 
the ground-state phase shift. 
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